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SMALL DEVIATIONS OF A GALTON-WATSON PROCESS WITH 

IMMIGRATION 


Nadia SidorovaQ 


Abstract: We consider a Galton-Watson process with immigration (Zn), with 
offspring probabilities (pi) and immigration probabilities {qi). In the case when 
Po = 0, Pi / 0, ( 7 o = 0 (that is, when essinf(Z„) grows linearly in n), we establish 
the asymptotics of the left tail P{VV < e}, as e 0, of the martingale limit W of the 
process (Zn)- Further, we consider the first generation 1C such that Zjc > essinf(Zx;) 
and study the asymptotic behaviour of 1C conditionally on {W < e}, as e 0. We 
find the scale at which 1C goes to infinity and describe the fluctuations of 1C around 
that scale. Finally, we compare the results with those for standard Galton-Watson 
processes. 


AMS Subject Classification: Primary 60J80. Secondary 60F10. 

Keywords: Galton-Watson processes, Galton-Watson trees, immigration, martingale limit, con¬ 
ditioning, lower tail, small value probabilities, large deviations. 


1. Introduction and main results 

1.1 Galton-Watson process without immigration 

Let {Zn : n > 0) be a supercritical Galton-Watson process with a non-degenerate offspring 
random variable X. We denote the offspring probabilities by {pk : k > 0) and assume that 
Po = 0. Further, we assume that the process starts with one ancestor, i.e. Zq = 1, and denote 
by a = EA > 1 the average offspring number. 

It is well known by the Kesten-Stigum theorem that under the condition EX log X < oo the 
martingale limit 

W= \im — 

n—>-oo qP' 

exists and is strictly positive almost surely. Moreover, the random variable W has a strictly 
positive continuous density, see [ 2 ] . However, only in a very limited number of examples can 
the distribution of W be computed explicitly, and one has to rely on asymptotic results to 
describe the behaviour of W. 

The left tail asymptotics F{W < e} as e 0 of IF has attracted a lot of mathematical atten¬ 
tion, both in its own right and in the broader context of small value problems, see [ 12 ]. Since 
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small values of the martingale limit W correspond to sub-average branching of the Galton- 
Watson tree, one naturally has to distinguish between the Schroder case when pi > 0 and the 
Bottcher case when pi = 0. In the Schroder case, small values of are achieved much more 
easily and, in particular, the minimal value of is equal to one with positive probability. 
In sharp contrast to that, in the Bottcher case the minimal tree grows exponentially. 

The Schroder case was hrst studied in [8], where it was shown that 

Ee“^'^ X s~^ as s —)■ oo, 

with r = — Shortly after that in [9] it was proved that the density w oi W decays at 
zero as w{e) x which easily implies 

P{1T < e} X e’" as £ I 0. 

The asymptotics for the density w was then rehned in [5] to w{e) ~ with an analytic 

multiplicatively periodic function L, which led to the corresponding improvement of the left 
tail asymptotics 

F{W < e} ~ L{e)e^ as e | 0 (1.1) 

with another analytic multiplicatively periodic function L. 


In the Bottcher case, the situation is more complicated as the tail of W decays exponentially. 
It was shown in [8] that the Laplace transform of W at inhnity has the logarithmic asymptotics 

logEe”^’^ X as s —oo, 

where (3 = with /i = min{n : Pn > 0} > 2 being the minimal offspring number. This 
suggested the logarithmic asymptotics 

logP{lT < £} X £“1^ as £ I 0 (1.2) 

and a more precise result was then obtained in [5j, namely, that 

logP{lT < £} ~ —M{e)e~^ as £ | 0, 

where M is an analytic positive multiplicatively periodic function. A numerical example with 
tiny but non-trivial oscillations in M was provided in [3], and an example with a constant 
M was given in m- Finally, the full left tail asymptotics was computed in [10] to be 

¥{W < e} ~ M(£)£ 2 (i-/j) exp { — } as e | 0, 

where both M and M are analytic, positive, and multiplicatively periodic. 

The precise form of the above asymptotics as well as the approach developed in [TO] made it 
possible to understand the influence of small values of W on the Galton-Watson tree. Let 


K = min{n : > pP"} 


be the hrst generation, where a vertex has more than the minimal number of offspring. We 
will call this event the hrst non-trivial branching of the tree. It was shown in [1] that, 
conditionally onW < e, the hrst branching time K will grow in the Schroder case as 


7s (e) 


log(l/g) 

log a 


(1.3) 





SMALL DEVIATIONS OF A GALTON-WATSON PROCESS WITH IMMIGRATION 


3 


and in the Bottcher case as 


7b (e) 


iQg(Vg) _ loglog(l/g) x 

log(a//i) log a ^ ^ ^ 


for some continuous multiplicatively periodic function H. However, more striking are the 
fluctuations of K. It was proved in [1] that in the Bottcher case there are no fluctuations 
at all, with K being equal to either [ 7 b(£)J or L 7 b(£)J + 1 with probability tending to 1, 
conditionally onW < e. This is no longer true for the Schroder case. Our hrst result below 
shows that the random variable K — 7 s(e) conditioned onW < e has exponentially decaying 
left and right tails. 


Theorem 1.1. In the Schroder case, as x ^ oo, 

liminf P{i^ — 7 s(£) > a: | IT < e} x limsupP{i^ — 7 s(£) > a: | IT < e} x 


liminf — 7s(e) < — a; | IT < e} x limsupP{i^ — 7s(£) < —x | IT < e} x 
=4-0 el-O 

where A = min{/c > 1 : p^ 7 ^ 0}. 

Remark 1. It was shown in [Tj that both tails of — 7s(£) conditioned on IT < e are not 
heavier than exponential but it was not known whether this estimate is sharp for either of 
the tails, and if so what the correspondent exponents and rates are. To hnd the actual tail 
behaviour, we had to control the error term of the asymptotics fll.ip . o 


Remark 2. Under the conditioning W < e there are two competing effects influencing the 
behaviour of K\ branching too early would lead to higher values of IT but on the other 
hand it would be probabilistically expensive to suppress branching for too long. Having hrst 
branching in generation roughly equal to 7 s (e) corresponds to the optimal trade-off between 
these two effects. The right tail of K — 'yfle) corresponds to a late branching, and its decay is 
given by the probability of having just one offspring in many generations, which is exponen¬ 
tial with exponent pi. The left tail corresponds to an early branching, which manifests itself 
in the appearance of extra (A — 1) offspring too early. The left tail is therefore controlled 
by the probability of keeping the sum of A (rather than one) i.i.d. copies of IT small, which 
explains the exponent p^“^ governing the left tail. o 


1.2 Galton—Watson process with immigration 

The remarkable difference in fluctuations of K in the Schroder and Bottcher cases is due to 
the fact that in the former setting the minimal tree does not grow at all, having just one 
offspring in every generation, whereas in the latter one the minimal tree grows exponentially. 
A natural question to ask is what happens if the process behaves similarly to the Galton- 
Watson process but its minimum grows linearly. Galton-Watson process with immigration 
is a natural example of such a process. 

Following the dehnition in [2], we £x a non-degenerate offspring random variable X with 
distribution (p^., fc > 0) as before, and an immigration random variable Y with distribution 
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(gfc : /c > 0). We define the Galton-Watson process with immigration {Z^ : n > 0) recnrsively 
by setting Zq = Yq and 

= X}"' + • • • + + W+i, n > 0, 

where all are independent and identically distributed with the same distribution as X, 
all Yj are independent and identically distributed with the same distribution as Y, and all 
X-"^ and Yj are independent. In other words, the Galton-Watson process with immigration 
{Zn) differs from the ordinary Galton-Watson process with offspring probabilities {pk) by the 
property that, in generation n, there is an immigration of a random number of individuals 
into the population governed by immigration probabilities {qk) and independent of the rest 
of the process. 

As before, we assume that po = 0. We also assume that pi > 0 as otherwise the linear effect 
of immigration will be negligible with respect to the exponential growth of the population. 
For the immigration probabilities, we assume that go = 0 in order to avoid the extinction 
and sub-linear growth of the minimal tree. 

We assume that 

EX log X < oo and E log Y < oo 

and denote a = EX, which is hnite by the hrst condition above and greater than one since 
Po = 0. It is a classical result, see [TB] for example, that under the conditions above the limit 

W = lim — 

n—>-oo qP 

exists and is positive almost surely. 

The following logarithmic left tail asymptotic was recently computed in [7]. As e 0, 

logP{yV’< e} ~ —(jlog^(l/£:), (1.4) 

where 

_ i^log(l/pi) 

(J — -^- 9 - 

2log^ a 

and 

u = min{i : g* > 0} > 1 

is the minimal number the immigration random variable can take with positive probability. 
As it was natural to expect, the left tail of W is thinner than that of the martingale limit W 
in the Schroder case fll.ll) but thicker than that of W in the Bottcher case fll.2l) . The above 
result was then generalised to multitype processes in [6]. 

1.3 Main results 

The aim of this paper is to hnd the full (non-logarithmic) left tail asymptotics of W at zero 
and to describe the time of the hrst non-trivial branching 

/C = min{?7, : Zn > + 1)} 

conditioned on W < e: in the limit as e J, 0. In particular, we want to compare the huctuations 
of /C around its typical growth with those for Galton-Watson processes without immigration 
in the Schroder and Bottcher cases. 
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Let u be the function defined implicitly in a right neighbourhood of zero by 

uj{e) — logCi;(e) + log log a = log(l/e). (1.5) 

In the sequel we will drop e in most notation and, in particular, in cu = a;(e), if there is no 
risk of confusion. We will also assume that e is sufficiently small. 

Theorem 1.2. yds e 0, 

P{W < e} ~ exp I — aoj'^ + uMi{uj) — - logo; + M 2 (a;)|, ( 1 . 6 ) 

where Mi and M 2 are bounded functions periodic with respect to uj with period logo. 

Remark 3. The leading term of the asymptotics fll. 6 l) agrees with the logarithmic asymp¬ 
totics oa since u{e) ~ log(l/£) according to fll.5p . o 

Remark 4. It is of course possible to express the asymptotics of P|W < ej in terms of £ 
rather than u but the formula would be too bulky and less transparent, and would require 
replacing uj by its asymptotic decomposition in e up to the hfth order term. o 


Now we turn our attention to the hrst branching time fC. Similarly to the immigration-free 
case, the second largest number of offspring 

A = min{i > 1 : pj > 0} 


will play an important role in describing the behaviour of /C. It turns out that, conditionally 
on W < e, the typical growth of /C is given by the scaling function 7 = 7 (e) dehned in a 
right neighbourhood of zero by 


7(£) 


lQg(l/g) I 
logo 


- 1 1 - 
-logo (A-1) log Pi- 


loglog(l/e). 


(1.7) 


The next theorem identihes 7 as the right scale and describes the fluctuations of /C around 
it conditionally on W < e. 


Theorem 1.3. There are constants ci,C 2 > 0 such that 

limsupP{/C > L7(e)J + x\yV < e] = exp { — cipf^^~^^^'\ (1.8) 

£4,0 

lira inf P{/C > L 7 (£)J -f a; | W < e} = exp { — C 2 pf^^~^'’^'\ (1.9) 

for all x e Z. 

Remark 5. Comparing 7 with 7 s in the Schroder case it is easy to see that the immigration 
can both force the process to start branching earlier and later. The former situation occurs 
if api~^ > 1 and the latter if ap^~^ < 1. This is not intuitively clear at the hrst glance as 
we would expect an earlier branching for small values of pi. However, the catch is that an 
early branching would require suppressing a larger number of the subtrees to ensure W < e, 
which is hard if pi is small. o 


Remark 6 . Theorem 11.31 describes the huctuations for every hnite a: G Z rather than just 
identifying the tail behaviour as x —?■ ± 00 . It immediately implies that the huctuations are of 
hnite order with a double-exponential right tail. This puts the Galton-Watson process with 
immigration strictly between the Bottcher case for the standard Galton-Watson process. 
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where there are no fluctuations at all, and the Schroder case, where the right tail decays 
exponentially by Theorem 11.11 At the same time, the left tail still has an exponential decay 
since 

1 f “(-^“1)^1 ^ —(A—l)a; 

1 - exp I - QPi I X Pi 

as X —?■ —oo for i = 1,2. Moreover, according to Theorem 11.11 the exponent of decay of the 
left tail is exactly the same as in the Schroder case, which shows that, unlike the right tail, 
the left tail is not affected by the immigration. o 

The comparison between the fluctuations in all three cases is given on the picture below. 



1.4 Ideas of the proofs 

It is well-known that the left tail of a positive random variable at zero is closely related to its 
Laplace transform at infinity. Following the approach suggested in [10], we use the precise 
inversion formula to obtain P{W < e} from the Laplace transform = Ee^^. This is a 

technically challenging step but, unlike the standard large deviations techniques, it is capable 
of providing the full asymptotics of P{W < e} rather than the logarithmic one. Then we 
use the immigration mechanism to relate with the Laplace transform = Ee^^ of 

the immigration-free Galton-Watson process. Further, we rely on the well-known Poincare 
functional equation 

( 1 . 10 ) 

where /„ is the n-th iterate of the generating function of the offspring random variable X, 
in order to understand the behaviour of (p{z) for large z through the asymptotics of /„ as n 
tends to infinity. 

The paper is organised as follows. In Section [2] we introduce notation and prove a couple of 
technical results for standard Galton-Watson processes. In Section [3] we establish the error 
term of the asymptotics fll.ip and use it to prove Theorem 11.11 In Section |3] we introduce 
notation relevant to immigration and get some preliminary asymptotic results for Galton- 
Watson processes with immigration. In Section [3] we prove two technical lemmas describing 
the behaviour of the Schroder function. In Section [3] we reduce the problem of describing 
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the behaviour of W and K, to that of understanding the left tails of a certain sequence of 
random variables (Vn), and establish some asymptotic properties of their Laplace transforms. 
In Section [7] we study the left tails of (V„) and, hnally, in Section [8] we combine everything 
and prove Theorems 11.21 and 11.31 


2. Notation and preliminaries in immigration-free case 


For any r > 0 , let Vr = {z E \ 1^1 < r} be the closed disc of radius r. Denote by 

OO 

= ZEVi, 

n=l 

the generating function of the offspring random variable X. Denote by 

ip{z) = z G C,7Zez > 0, 

the Laplace transform of the martingale limit W of the corresponding Galton-Watson tree. 


Let fo{z) = z and, for each n > 1, denote fniz) = f{fn-i{z)). Since the behaviour of the 
iterates of f{z) for large n is mainly determined by its leading term piz, it is convenient to 
use the decomposition 

n—1 

fn{z) = f{fn-l{z)) =Plfn-l{z)(l+p]:^J2Plfn-\i^)) =Pi^YIM^)^ ( 2 - 1 ) 

1>1 j=0 

where the functions Aj are dehned on "Di by 

A,{z) = l+p:['J2p,f‘p{z). (2.2) 

1>1 


For each \z\ < 1, denote 


S{z) = zYlAj{z). (2.3) 

j=0 

It is well-known (see [H Lemma 3.7.2 and Corollary 3.7.3]) that this inhnite product converges 
uniformly on each r G (0,1), and the function S is called the Schroder function. It is 
easy to see that 

= lini Mf). (2.4) 

In particular, on each r G (0,1), the Schroder function S is bounded and Aj(z) —?■ 1 
uniformly as j —?■ oo. Denote 


B = {zEC:\z\ <1,F(2)^0}. 
For each n > 0 and \z\ < 1, denote 

Rn{z)=S{z)-p^-Uz). 


(2.5) 


The asymptotic behaviour of Rn will be crucial for our analysis. In the remaining part of 
this section we prove some elementary properties of the functions S, Rn, and An. 


Lemma 2.1. faj | 5 '( 2 ;)| < Fd^^l) for any \z\ < 1. 

(b) The functions s i-A S{s) and s i-A S{s)/s are increasing on [0,1). 
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Proof, (a) This follows from the same property for each Aj which, in turn, follows from the 
same property for each fj. 

(b) Obviously, it suffices to prove the second statement only. Observe that S'(s)/s is a product 
of Aj{s), and each Aj is increasing since each fj is increasing. □ 

Lemma 2.2. (a) Let r G (0,1). Then there are constants ci,C 2 > 0 such that 

An{z) - 1 ~ 

and 

Rn{z) ~ C2S‘^{z)pf^^~^^ 


as n ^ oo uniformly on Vj. A B. 

(h) Rn{z) =0 if z E Vr\B for all sufficiently large n. 

Proof, (a) Every convergence and equivalence mentioned in the proof below is meant to be 
uniform on Pr A B. 

Using fl2.ip and fl2.3p we obtain 


Rn{z) = Pi ( JJ Ak{z) - 1 j. 

k=n 

Since Aj(z) ^ 1 as j —?■ oo and using fl2.4l) we have 


OO 


Rn{z) ~ S{z) 


exp 


J2^ogAj{z) 

j=n 


1 


( 2 . 6 ) 


Observe that Aj{z) ^ 1 eventually on B. Indeed, it follows from (12. 4 p that fn{z) —)■ 0 as 
n ^ OO. Hence the hrst term of the sum in (12.2p dominates over the remaining terms 


i>\ 




oi\fYi^)\). 


(2.7) 


It remains to notice that fn{z) ^ 0 eventually by fl2.4p and use fl2.2p together with fl2.7p . 


Now the first statement follows from 

log Aj{z) ~ Aj{z) - 1 ~ ~ PxPi^~^^~^S{z) 

with Cl = p\/Pi, where the middle equivalence is implied by fl2.7p and the last one by fl2.4p . 
Now we have 


^logHj(z) ~ ~C2^(z)pi^^ 

j=n j=n 

with some C 2 > 0. Substituting this into (12.61) we obtains the required asymptotics. 

(b) It is easy to see from fl2.3p that S'(z) = 0 implies z = 0 or Aj{z) = 0 for some j. Then 
fn{z) = 0 eventually by (12.1 p and so Rn{z) = 0 as well. 

□ 
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3. Fluctuations in the Schroder case 


In this section we prove Theorem 11.11 The proof for the right tail of iF — 7 s will be rather 
straightforward. The left tail, however, is controlled by the second term of the asymp¬ 
totics fll.ll) . which we estimate in the lemma below. 

Lemma 3.1. In the Schroder case, 

L{e)e^ - ¥{W < e} x 

as e 0, where r = —7^ and L is an analytic multiplicatively periodic function on (0, cxo) 
with period a. 

Proof. Since W is almost snrely positive and has no atoms, its left tail can be compnted by 
the inversion formula 

1 /‘°° 1 - e-"^ 


P{1F <e} = 


2tt 


IT 


■ ip{—iT) dr. 


(3.1) 


Recall the dehnition fll.3p of 7s. Changing the integration contour from the vertical axis to 
the vertical line passing through and substituting r = we obtain 

^00 ( 1 —it) ^ 


1 r 

F{W <e} = - f 


1 — it 


(p((l — dt. 


(3.2) 


By the dehnition of 7 s we have = a Further, the Poincare functional equa¬ 

tion fll.lUp and fl2.5p imply 

9 ?(( 1(<^(1 - R)) = - R)) - - *^))- ( 3 - 3 ) 

Substituting this into fl3.2p and taking into account that = e'^Pi by the dehnition 
of r, we obtain 


W <e} =e' 


Pi «i) _ 


27r 


1 — it 


-R)) -Rl7sJ (</?(! dt- (3-4) 


This representation naturally splits the the left tail probability into the leading term corre¬ 
sponding to S and the error term corresponding to R[ 77 . Namely, we dehne 


L{e) = 


Pi 


~{7s} /*oo a it) -1 

1 / t; 1 


27T 


S{(p{l — it)) dt. 


(3.5) 


1 - it 

. -00 -*■ 

In order to prove that the integral is hnite we hrst observe that the ratio under the integral 
is bounded. Then we get by Lemma 12.11 


|S(v>(l - !i))| < S(|»)(l - it)\) < -^^^lv>(l - *0 


(3.6) 


and use m Lemma 16] which claims that |<p(l — it)\ is integrable over M with respect to t. 
Hence the function L is well-dehned, bounded, and multiplicatively periodic with period a 
since 7s(£) — 7s (as) = 1 G Z. In particular, once we have shown that the error term is 
negligible, it will imply that the function L must be the same as in fll.ip and hence real¬ 
valued and is analytic. 

For the error term, we use Lemma [2.21 with r = (p(l) to get 

foo a) _ 


p-{7T /■oo ga-{T'H(l-ip _ ^ 


27r 


1 — it 


^L 7 sJ( 9?(1 -it))dt^e 


tX 


it 


S {if {1 — it)) dt, 
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and it remains to show that the integral on the right-hand side is bounded away from zero 
and inhnity. 

To do so, consider two independent random variables Wi and W 2 with the same distribution 
as W. Observe that 


X P{lTi < e/2, W 2 < e/2} < ¥{Wi + W 2 <e} < ¥{Wi <e,W 2 <e}- (3.7) 

Similarly to fl3.ip . (13■2p . (13■3p . and (13.4^ we have 

1 /■°° 1 - 


F{Wi + W2 <£} = 


2vr 


ZT 


■ (p {—ir) dr 


'■00 go 

I- it 


[^( 99(1 - it)) - - it))Y dt. (3.8) 


By Lemma 12.21 and boundedness of S we have 

roo (l—it) _ ^ 


1 - 

r*oo (1—'it) ^ 


- it))Rpy^^{p{l - it)) dt = b) / \s^{p{l - it))\dt, 

J —00 

POO 

=0(e2^(^-b) / |,S^((y9(l-it))| dt. 


/_oo ^-it . _ 

Since S is bounded, both integrals on the right-hand side are hnite by the same argument as 


above in (13.6p combined with [Tnl Lemma 16]. This implies that both terms on the left-hand 
side are o(l) and hence the main term of the asymptotics (13.81) is given by 

coo it) _ 


¥{Wi + W 2 <e} X£ 2 r 


1 — it 


o ((p(l — it)) dt. 


Now fl3.7p implies that the above integral is bounded away from zero and inhnity. 


□ 


Proof of Theorem li.il It suffuses to prove the theorem for a; G Z. Observe that the condition 
K > 7 g -|- a: is equivalent to having just one offspring in the generation |_ 7 sJ -|- x. On this 
event, the condition W < e is equivalent to IT < where W is the martingale limit of 

the Galton-Watson subtree generated by that offspring. Hence W has the same distribution 
as W but is also independent of the event = 1. Using fll.3p we obtain 

F{K > 7g + X, lU < e} = P{^L7 bJ+^ = 1, W < = pL7sJ+xp|^ ^ j. 


Combining this with the left tail asymptotics fll.ip and using pf = e'^ we get 


[ 7 aJ +X 

P{iL > 7s -F X I lU < e} ~ P{lU < 

J-j\S 


,-{7s}+a:j _ /El 


-{7s}+x 


L(a-UU) 


F{W < 


-{7s}+a:' 


as £ j, 0 since L is multiplicatively periodic with period a and e = 0 “'’'=. Observe that, for a 
hxed X, the expression on the right-hand side only depends on {7s} and so is multiplicatively 
periodic in e. This implies that 


lim inf P|iL > 7 s -|- x | fU < 
limsup P{iL > 7s -I- X I lU < e} 

eiO 


Pi 

^ L(a-“i) ^ ^ 

—Q ;2 

Pi 

^ L(a-"2) ^ J’ 
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for some ai,a 2 G [0,1]. It suffices now to show that uniformly for all a G [0,1] 


Pi 




and 


as X ^ oo. The hrst identity easily follows from the fact that 

¥{W < > F{W < a-^+^} ^ 1 


as X —)■ cxo. For the second one we observe that 0 as x ^ cxo and so we can derive 

the above asymptotics from the left tail asymptotics of the martingale limit obtained in 
Lemma 13.11 As L is multiplicatively periodic with period a and x is an integer we have 




L(a-«) 


} X - P{1F < a 






(A—l)a: 
Pi 


since a '^ = pi. 


□ 


4. Immigration 

In this section we introduce notation relevant to the immigration and prove some preliminary 
results which will be necessary to deal with it. Denote by 

CX) 

h{^) = '^(lnz'^, zeVi 

n=u 

the generating function of the random variable Y. We will see in Section [6] that in order to 
hnd the left tail asymptotics of W we will have to control the products 

N 

IlM/nW) (4.1) 

n=\ 

for G N. Since fn{z) will typically tend to zero as n —)■ cx), the function h{z) will essentially 
behave according to its leading term QuZ'^. This observation suggests using the decomposition 

hUn{z)) = q^Z{z){l + q~^'^qif~''{z)^ = qJ^{z)Bn{z), (4.2) 

1>U 

where the functions Bn are dehned by 

= l + zEV,, (4.3) 

1>U 

Combining fl2.ip and fl4.2p we obtain 

n—1 

h(Uz)) = n A(^)- 

j=0 

For any r G (0,1) and 9 G (0, vr/2), denote 

'^r,e = {z E Vr : z ^ 0,\ aTgz\ < 9}. 


(4.4) 
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Lemma 4.1. Let r G (0,1). As n ^ oo, 

n 

X{b,{z)^B{z) 

i=i 

uniformly on Vr, where B is a bounded holomorphic function. 

Further, there exists 6 G (0,7r/2) such that on Vrp 

(a) S{z) ^ 0 and, in particular, A^^z) ^ 0 for all k; 

(b) B(z) ^ 0; 

i=i 

uniformly as n ^ oo, where C is a hounded holomorphic function. 

Proof. Using fl4.3p and monotonicity of /„ as well as that it is bounded by 1 on (0,1) we get 
\Bj{z) - 1| < qf^^qi\fj~’'{z)\ < qf^fj{r) ~ qf^S{r)p{ (4.5) 

for all z ^ Fr uniformly. Since the sum over j of the expressions on the right hand side is 

n 

hnite, Bj{z) converges uniformly on Dr, and hence B is holomorphic and bounded. 
i=i 

Let us now prove the second part of the lemma. 

(a) Observe that S is holomorphic on Dr as a uniform limit of holomorphic functions. Hence 
it can only have isolated zeroes, and it suffices to observe that S'(s) 7 ^ 0 for s G (0,1) and 
choose 6 sufficiently small. 

(b) Similarly, B can only have isolated zeroes, and it is easy to see from 04.31) that B{s) > 1 
for s G [0,1]. Hence one can choose 6 small enough so that B has no zeroes on Drp. 

(c) First we observe that by (a) the product is well dehned for a sufficiently small 6. Second, 
by Lemma [2.21 we know that 

j\ogAj{z) ~ cS{z)jp]^^~^^ 

uniformly as j 00 , where S{z) 7 ^ 0 by (a) if 6 is small enough. Again, the sum over j of 
the expressions on the right hand side is hnite, and we can apply the same arguments as for 
the hrst product. □ 


Lemma 4.2. (a) Let r G (0, 1). There is c> 0 such that, for all s G (0,r], N and z G Dg, 

N , , 




n=l 


7 iziV(iV+l) 


(4.6) 


(b) Let r G (0,1) and let 9 be chosen according to Lemma f.l. Then 

N 


W^h(U{z)) ~ F{z)q^p, " /yfi(7 


(4.7) 


n=l 

as N ^ 00 uniformly on Drp, where F is a bounded function on Drp, which is nowhere equal 
to zero. 
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Proof. It follows from fl4.4p that, for all iV > 1, we have 


N 


Ylh{fn{z)) = 


uN{N-\-l) 
2 


N n—1 


N 


n=l 


nnAw)(n«” 

n=lj=0 


n=l 


(4.8) 


(a) Observe that the last term is uniformly (in s and z) bounded by Lemma 14.11 Further, it 
follows from fl2.2p that < Aj{s) for all j on Vg. Hence 


N n—1 


N n—1 


N 




'■■nriAw <iT‘'"-‘nnAw 

n=l j=0 n=lj=0 


= \z\s' 


vN- 


i=o 


pN-j) 




where we included j = N into the product since the corresponding term equals one. As all 
Aj{s) > 1 we can drop j in the exponent which together with fl2.ip implies 


..uN 


N n—1 

nriT 

n=lj=0 




(b) It follows from Lemma 1401 that all Aj{z) ^ 0 and so we can rearrange the middle product 
in 04.81) and use 02.ip to get 


N n-1 N N 

n n n 




n=l j=0 j=0 j=0 

Now the statement of the lemma follows from Lemma [4.11 with F{z) = B{z)C‘'{z). □ 


5. Finer properties of the Schroder function S 


Lemma 5.1. The function s i—)■ logS'(s) is well-defined and analytic on (0,1), and 


(log^(s))' > 1/s 


for all s G (0,1). 


Proof. Since S'(s) > s by 02.3p the logarithm is well dehned in a complex neighbourhood of s, 
and the analyticity follows from log S being holomorphic there. Differentiating the uniform 
limit of holomorphic functions 02.4p we obtain 


(log^(s))' = 


S{s) 


= lim 

n—^oo 


f'nis) 

fn{s) 


> 1/s 


since this inequality is true term by term for ff and /„. 


□ 


Denote 


Lemma 5.2. 


i/(s) = logF( 9 ?(s)), s>0. 

The function -0 is analytic, '4’'\s) > 0 for all s, and 

lim ijifs) = 0 and hm'^'(s) = —cxo. 
s^oo s4-0 


(5,1) 
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Proof. The function if is analytic as a composition of analytic functions. Compute 


lf'{s) = 


r{s) = 


S{ip{s)) 

SHpis)) 


It was shown in [TUI (75)] that 9 ?"(s) 9 ?(s) > (</?'(s))^ for all s > 0. Further, S is positive 
according to (12.31) and S' is positive by Lemma 15.11 This implies 


r{s) > 


{S'ipis)) + <fi(s)S"((fi(s)))(ip'(s))^S(ip(s)) - ip(s)(fi'(s)(S'(ip(s))y 


S^(cp(s))cp(s) 
Using [ini (63)] for /„ instead of / we obtain 

'^fn(s) , ^ Q 


(5.2) 


fn{s) 

for all s G (0,1) and all n, which extends to 

'sF'(s) 


> 0 


for all s G (0,1) since one can differentiate uniformly convergent series of analytic functions. 
This implies 

(F'(s) + sF"(s))^(s) > s(F'(s)) 2 
for all s G (0,1). Substituting this into fl5.2p we get 

F2(93(s)) 

for all s > 0 since (p'{s) < 0. 

To prove fl5.ip we observe that by the Poincare functional equation and fl2.4p we have 
V'(sa") = \ogS{fn{(pis))) = \og{pfS{ip{s))) = if{s)+n\ogpi 

and so 

a^if'^sa") = if'is) 

for all s > 0 and all n. Since if" is positive if' is decreasing and so 

lim if'{s) = lim if'{a") = if'{l) lim aP" = 0, 

s—^oo n—>-cx) n—^oo 

hm'0'(s) = lim if'(a~") = if'{l) lim a" = —oo, 

s4,0 n—^oo n—^oo 


as required. 


□ 
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6. The random variables and their Laplace transforms 


Denote by 


= Ee 2 ; G <C,7lez > 0. 


( 6 . 1 ) 


the Laplace transform of the random variable W. Observe that the m-th generation of the 
Galton-Watson tree with immigration can be written as 

Vo 

= ( 6 . 2 ) 

i=\ 

where where is the Galton-Watson process corresponding to the Lth immigrant in the 
generation zero, and Z is the Galton-Watson process with immigration starting with the 
immigrants of generation one. It is easy to see that the process Z and all processes are 
independent, and Z has the same distribution as Z. Dividing by aZ and passing to the limit 
we obtain 

Vo 

W = ^Wi + a-^W, (6.3) 

i=l 

where W are the martingale limits of the processes and VV is the limit random variable 
corresponding to Z. Glearly, W and all Wi are independent and have the same distribution 
as W and W, respectively. This implies 

ip*{z) = h{(p{z))(p^{za~^). (6.4) 

Further, for any k > —1, denote by 

u{k+l) 

Vfc = a-^ (6.5) 

i=l 

where W-^\ 1 < i < + 1), are the martingale limits of independent Galton-Watson 

processes indexed by the hrst v immigrants of all generations between 0 and k, and 
is the limit random variable of the independent Galton-Watson process with immigration 
starting with the immigrants of generations strictly after k. Similarly to fl6.3p we observe 
that 


W = Vfc on the event {/C > k} 

since the minimal number of individuals in generation k is ^{k + 1) given by u immigrants in 
each generation and having just one offspring each. Denote the Laplace transform of Vk by 

(Pk{z) = 2 ; G C,Tlez > 0. 

It follows from fl6.5p that 

ipk{z) = ( 6 . 6 ) 

Observe that V_i = W and = ip^. We also denote Z_i = 0. 

For each fc > — 1 we have 

f** fc (k I 1) 

P{/C > k,W <e} =F{Zk = u{k + l),Vk<ej= < e}. (6.7) 

This means that our main aim now is to understand the left tail probabilities of Vk- Those 
corresponding to fc = — 1 will give us the left tail asymptotics of W, and those with k = 
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7 (e) + X will control the fluctuations. Recall that the left tail of Vk is closely related to the 
behaviour of the Laplace transform ^pk for large values of the argument. The next lemma 
enables us to understand it through the asymptotic properties of the iterations /„ and the 
branching mechanism h in the same spirit as the Poincare functional equation does it for a 
standard Galton-Watson tree. 


For any k and N, denote 

Further, for any z E B, denote 


i/JV{JV+l)-i/fc(fc+l) 

— Qu Pi 




:,n{z) - (^1 


RN-k{z) \ 

~W) 


Lemma 6.1. (a) Let r G (0,1). There is c> 0 such that 

\M^a^)\<c^-^Ck,NS{sr^. 

s 

for all k < N, s E ( 0 ,r], and all z satisfying \p{z)\ < s. 

(b) Let r E (0,1) and let 6 be chosen according to Lemma \4.1\ Then 

ifkiza^) ~ ^k,N{T{z))p>^{z)F{p>{z))Ck,NS{ip{z)y^ 
as N — k ^ 00 uniformly on {z : ip{z) E Vrfi}- 
Proof. Using (16.bh and iterating it according to (16.dh we obtain 

N-k-l 

(pk{za^) = ipyz)ip{za^~y’'^^^^'> h[ip{zay). 

The Poincare functional equation fll.lOp implies 

ifkiza^) = (pyz)ffl>^'l:^\(p{z)) JJ h{fn{p{z))). 


n=l 


N-k-l 


n=l 


(a) Estimating the Laplace transform by one, the second term of fl6.10p by 

\fN-kiTiz))\ < fN-ki\Tiz)\) < fN-kis), 

and the product by Lemma 14.21 we have 

\(Pkiza^)\ < c- 


Mz)\ . 


Pi 


fN-kis) 


( 6 . 8 ) 


(6.9) 


( 6 . 10 ) 


with some c > 0. Taking into account 

fN-k{s)<P^-^S{s), 

which follows from fl2.1l) and fl2.3p we obtain the required bound. 

(b) Using the asymptotics for the product in (I6.10p obtained in Lemma 02] we get 

ifkiza^) ~ T*{z)F{ip{z))qy-’'-^p^ ^ fN-kipiz)) 

as — fc —)■ 00 uniformly on {z : (p{z) E T>r,e}- Taking into account 

fN-kipiz)) = Pi~’'S{(p{z)) -p^-’^RN-kipiz)) 
and observing that S{(p{z)) y 0 on {z : (p{z) E T’r.e} by Lemma ITT] we arrive at 06.91) . □ 
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7. Left tail of Vk 


The aim of this section is to compute the left tail asymptotics for Vk for two types of k. The 
hrst case is simply 

k = -1. (7.1) 

Combined with fl6.7p . this would give us the left tail asymptotics of W and prove Theorem 1 1.21 
The second case is 


k = k{e,x) = \_^{e)\ + X (7.2) 

for a hxed integer x. This is needed to prove Theorem 11.31 It turns out that both cases can 
be handled simultaneously so in this section we always assume that k satishes either fl7.ip 
or (17^ . 

Let p be the function dehned implicitly in a right neighbourhood of zero by 

= e. (7.3) 

It is easy to see that 


u = ploga. 

and that the hrst three leading terms of the asymptotics of p are given by 


P(^) = 


log a 


log(l/£) + loglog(l/e) - (1 + o(l)) log log a 


(7.4) 


(7,5) 


It is worth mentioning that the dehnition fll.7p of 7 manifested itself from the condition 

X 1 as e 10, (7.6) 

which will prove to be crucial later on. 

Now hx 


N = N(e)=[p(e)\. 

and choose u = u{e) in such a way that 

ixijj'iu) = (7.7) 

This is possible by Lemma 15.21 since ip' takes all negative values. Moreover, since the right 
hand side of fl7.7p is bounded between —1 and —1/a and ip' is decreasing by Lemma 
there exist positive and independent of e constants and u* such that u G for all 

sufficiently small e. 


Similarly to fl3.1l) . the lower tail of Vk can be computed by the inversion formula 

1 1 - 


l?{Vk<e} = 


271 


IT 


</?fc(-ir) dr. 


Similarly to fl3.2p . we will move the integration contour far to the right and rescale the 
integration accordingly. However, we have to do it more carefully. Namely, we replace the 
vertical coordinate axis by the vertical line passing through ua^ and substitute r = ta^. 
This gives 




^00 (u—it) 


Vk{{u — it)a^) dt. 


u — it 


(7.8) 
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In order to understand the impact of the term in the integral, we will combine Lemma l 6 TT] 
with the saddle point approximation for S. Choose r in such a way that < r < 1 

and choose 6 according to Lemma [4.11 The composition S o ip maps the interval [u^,,u*] to 
[S'((p(m*)), S{(p{u^,))] since p is decreasing on [0, oo) and S is increasing on [0,1) by Lemma [2Tl 
Hence we can choose Pi > 0 small enough so that p maps [u^,, u*] x [—/3i, /9i] to Vrp and Sop 
maps X [—/9i,/3i] to {z G C : TZez > 0}. Now the function log(S' o p) is well-dehned 

on X [—/3i,/9i] and its third order derivatives are bounded. Expanding into the Taylor 

series 

log S{p{u - it)) = logS'((p(M)) - itp'iu) - —i>"{u) + 0{t^) (7.9) 

as t —?■ 0 uniformly in e and using the fact that p" is positive by Lemma 15.21 we choose 
P G (0, Pi) in such a way that 

ne[\ogS{p{u - it))] < log S {p{u))-^—i>”{u) (7.10) 

for all t G [—P, P] and all e small enough. Finally, we choose a = a{e) so that a 0, 

0 and a^N/logN ^ oo. (7-11) 

We assume that e is small enough so that a < p. In Lemmas 17.2117.31 and 17.41 below we will 
compute the main part of the integral fl7.8l) coming from integrating over [—a, a], and show 
that the integrals over the remaining parts {|t| G [a,/3]} and {|f| > /5} are negligible. 

However, before turning our attention to the integral fl7.8p we compute the asymptotics of 
'^k,N, which plays a crucial role for pk{{u — it)a^) according to Lemma [ 6 Tl 


Lemma 7.1. Let r G (0,1) and let 6 he chosen according to Lemma \pPl 
(a) T„i^Ar(z) ~ 1 as e 0 uniformly on T>rfi. 

(h) Let k be of the form fl7.2l) . Then there is w = w{e) bounded away from zero and infinity 
such that 

^k,N{z) ~ exp { - S{z)] 

as £ 0 uniformly on Vrp . 


Proof. Comparing the dehnition of 7 in fll.71) and the asymptotics fl7.5p of p it is easy to see 
that N — k^oonselO in both cases. Observe that S{z) 7 ^ 0 on Vrp by Lemma 14.11 
Lemma [2.21 now implies that 


RN-k{z) 


Cl S{z)pt^'>^^-^^ 


with some Ci > 0 and hence 


logTfc^TV- CiuNS{z)p[^ 


as £ 4 , 0 uniformly on 

(a) If /c = — 1 then we use N ~ to obtain 

= 0 . 

£ 4.0 

(b) For A; = 7 — { 7 } + x we get 
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Since ~ C 2 according to fll.7p and the three leading terms of p given by fl7.5p . we 

obtain the required asymptotics with w = w{e) = ciC 2 P^i □ 

Lemma 7.2. Suppose that k is of the form fl7.1l) or (17.2|) . Then 

1 (u—it) 


2tt 


u — it 


Tk{{u — it)a^^ dt 


^k,NCk,N GW + uN log S{p{u)) - ^ log ivj 


as £ J, 0, where 


^k,N = ^k,N{^) — - l====^k,N{T{'^))- 


U 




\u) 


Proof. By our choice of (3 we have p{u — it) G Vrp for all t G [—/3,/3]. By Lemma l6TT] we 
obtain 

(Pk{{u - it)a^) ~ ^k,N{T{u - it))p^{u - it)F{p{u - it))Ck,NS{p{u - it))"^ (7.12) 

as £ J, 0 uniformly for all t G [—/9, /?]. Taking into account that a 0 we have 
ifk{{u - it)a^) ~ ^k,N{T{u))p^{u)F{p{u))Ck,NS{p{u - it))''^ 
as £ 0 for all t G [—a, a] since '^k,N is regular enough by Lemma ITTI We obtain 


r*a {u—it) 

U — it 


Tk{{u — it)a^) dt 


Ck,N'^k,N{T{,u)) 


p,{u)F{p{u)) r N 


u 


^^-^^^S{p{u-it)y^ dt. (7.13) 


For the integral above, we use the Taylor expansion (17.9p to obtain 

f e^^''^^-^^^S{p{u-it))''^ dt 
J —OL 

~ exp + i/A^logS'((^(M))} f exp I — iti^ea^ + uNy^u)) 


FvNif'iu) 


dt, 


where we have also used ae^“^ W 0 from fl7.1ip to get rid of the negligible terms in the 
Taylor expansion. Observing that by fl7.3l) and (17.71) 

ea^ + vNif'iu) = + vpfj'iu) — v{p}y{u) 

= + ny'iu)) — z/{p}'0'(m) = —v{p}y'{u) 


we obtain 




exp [uea^ + uN log S{(p{u))} / exp <itiy{p}ifj'{u) 


u - 


FvNy'iu) 


dt. 


(7.14) 
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Substituting r = we get 


exp 


[u] 


dt 




vN'lp"{u) J -a^JvNiP"{u) ^ \/vN'lp''{u) 2 J ^vNip"{u) 

since, as £ J, 0, the interval of integration increases to M by fIT.lip and the first order term 
tends to zero (see for example [U Lemma 11]). Combining this with fl7.14p we obtain 


dr 






exp + z/iVlogS'(v9(M)) — ^ logiV [>. 




u 


Together with fl7.13p this proves the required asymptotics. 
Lemma 7.3. Suppose that k is of the form CHi or (17.21) . Then 


□ 


'\t\&[a,l3] 


(u—it) 

u — it 


■ Tk{{u — it)a^^ dt 


< 


Cfc^Arexp + uN log S{(p{u)) — logA^j, 


for all e small enough. 


Proof. Observe that (I7.12p is in particular true for all t such that |t| G [a,/S]. Since \u — 
it\ > M*, |(p*(m — it)| < 1, F is bounded by Lemma 1121 and '^k,N is uniformly bounded by 
Lemma 17.11 there is a positive constant c such that 

(u—it) 

■ Tk{{u — it)a^) dt 




u — it 


< cCkM 


(u—i 


0|e[a,/3] 


‘>S(tp(u-it))''^ dt 


(7,15) 


For the integral above, we use the Taylor bound (I7.10p to obtain 


J\t\e[a,l3] 

< exp [uea^ + uN log S{(p{u))] 


< 2/9 exp + uNlog S{(p{u)) — 


' exp 
|t|e[a,/3] 

a^vNil)"{u) 


t^vNif'iu) 


dt 


(7,16) 


(7.17) 


Taking into account the fact that a^N/ logiV —)■ oo according to (17.lip and combining (I7.15p 
with fl7.17p we obtain the desired estimate. □ 

Lemma 7.4. Suppose that k is of the form dmi or fl7.2p . There exists 5 > 0 such that 


Tl>/3 


.^ea^ {u—it) 
U — it 


Tk{{u — it)a^) dt 


< C'fc,Aexp uea^ + z/iVlog S'(v9(m)) — 5uN 


(7.18) 
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and 


— it 


— ^Pk{{u — it)a^) dt 


< Ck,N exp < vN \og S{(p{u)) — duN 


for all e small enough. 


(7.19) 


Proof. Similarly to the proof of m Lemma 16], we use the fact that, for each v G [u^,u*], 
t i-A (p{v — it)/ip{v) is the characteristic function of some absolutely continuous law (Cramer 
transform), the continuity of the mapping {v,t) eA (p{v — it)/(p{v), and the compactness of 
[n*,n*] to conclude that there is a constant 7] such that 

\g:>{u — it)\ < (1 — ri)i.p{u) for all |f| > f3 (7.20) 

and all e small enough. Using Lemma 16.11 with s = s{e) = (1 — ri)(p{u) < (p{u^) < r and 
taking into account that s > (1 — rj)^p{u*) and so is separated from zero, we obtain, with 
some positive constant Ci, 

\ipk{,{u - it)a’^)\ < ci\ip{u - it)\Ck,NS{{l - vi)Lp{u)y^. 

By Lemma 15.11 and the mean value theorem 

log ^((1 - ri)g){u)) < log S((f(u)) - t], 

which implies 

\(Pk{iu - it)a^)\ < ci\g){u - it)\Ck,NS{g){u)y^. 

Substituting this estimate into the integral flT.lSp and using \u — it\ > u* we obtain 

/ - r— Pk{{u - it)a^) dt 

J\t\>g u-it 

< C 2 Cfc,AT exp log S'(v 9 (m)) — 771 /I f \ip{u — it)\dt (7-21) 

^ J —00 

with some C 2 > 0. It was shown in m Lemma 16] that the integral above is uniformly 
bounded. This implies the required estimate with some d < rj. 

The estimate fl7.19p is obtained in the same way as fl7.18p with the only difference that the 
terms (u-it) uea^ are omitted. □ 


8. Proofs of the main theorems 

In this section we establish the joint probability fl6.7p and then prove Theorems 11.21 and 11.31 
Proposition 8.1. Suppose that k is of the form 07.1 p or 07.2p . ^Is £ j, 0, 

F{]C > k,yV < e} q^p^ ^ ^k,NGW + oN log S{ip{u)) —-log N j. (8.1) 
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Proof. Combining fl6.7p and (17.81) we obtain 

Vk(k+1) 1 /■°o ea^{u-it) _ i 

P{/C > fc, >V < £} = ^ A-/ -— T - ( Pk{{u - it)a^) dt. 

ZTT J _r^ XL %L 


Let us split the integral into the sum of the three integrals corresponding to keeping 
in the numerator and integrating over [—a, a], {t : |t| G [a, /?]}, and {t : |t| > /?}, respectively, 
and the integral corresponding to keeping —1 in the numerator and integrating over M. 
Lemma 17.21 gives the asymptotics of the hrst integral, while Lemmas 17.31 and 17.41 imply that 
the remaining three integrals are negligible since is bounded away from zero and inhnity 
by Lemma 17.11 Substituting the asymptotics for the hrst integral given by Lemma 17.21 we 
arrive at the required formula. □ 


Proof of Theorem \1.2[ We use Proposition 18.11 with /c = — 1 as 

P{>V <e} = P{/C > -1, >V < e}. 

By Lemma [7. II we have 

P*{u)F{p{u)) 




N 


Uy^2Tri'fj'‘ 


[U 


Using (17.3p we get 


uea^ = uea^ = upa 


Further, 


uN log S{ip{u)) = uplogS{ip{u)) - u{p} log S{p>{u)) 
logiV = logp + o(l) 

du = exp {p log - {p] log ] 


uN{N+l) 

Pi ^ = exp 


u 


l^(p^ + Pl-2{p})+{pp-M)}, 


Substituting all of the above into (18.11) we obtain the required asymptotics (11.61) with 

Mi{u) = - - Mi( ^ 


log a V log a / ’ 

M2(a;) = i log logo+ M2fr~ /’ 

2 V log a / 


where 


Mi{p) = _ 2{p}) +ua + u log S{(p{u)) + log q^ 


M 2 {p) = log 


P*{u)F{p{u)) 


+ - {p}) -Kp}log>^(‘d(w)) - {p}^ogq^, 


which are bounded periodic functions of p with period one since u is a bounded function of 

{p} by ([72D- □ 
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Proof of Theorem 2. Let x G Z. Using Proposition 18.11 both with k = —1 and k = 
[ 7 J + X we obtain 


P{/C > L 7 J + X I W < e} = 


P{/C > L7J + X, W < e} <hL^j+ 3 ;,Ar 
P{/C > -1,>V < e} 


by Lemma [ 7 .1[ This implies (II.Sh and (ll.Qj) with 
Cl = i^lirninf and 


— (A—l)a; 

exp < — wup-^ 


S{p{u)) 


C 2 = i^limsup 

el-O 


k{e) = 


which are both positive and hnite since S o ip is continuous on [m*, u*] and w is bounded away 
from zero and inhnity by Lemma I7.1I □ 
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